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We investigate the liquid-vapor interface of the restricted primitive model for an ionic fluid using a density 
functional approach. The applied theory includes the electrostatic contribution to the free energy functional 
arising from the bulk energy equation of state and the mean spherical approximation for a restricted primitive 
model, as well as the associative contribution, due to the formation of pairs of ions. We compare the density 
profiles and the values of the surface tension with previous theoretical approaches. 
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1. Introduction 

According to the simplest model of ionic fluids, called the restricted primitive model (RPM), 
the fluid consists of charged hard spheres (at total density p) of equal diameter, d, half of which 
carry a charge of +q, while the other half carry a charge of —q. These spheres are immersed in a 
dielectric medium of dielectric constant e. Reduced thermodynamic quantities, temperature and 
density, appropriate to this model are defined as follows T* = kTed/q 2 and p* = pd 3 . 

Perhaps, the first successful statistical-mechanical description of this model was proposed by 
Debye and Hiickel [l[ about ninety years ago. The Debye-Hiickel approach predicted the existence 
of the first-order transition for RPM and was subsequently applied in several works, e.g., to describe 
phase behavior of three-component ionic fluids [2j]. More recently, the bulk RPM was studied by 
means of the Ornstein-Zernike (OZ) equation supplemented by the mean spherical approximation 
(MSA) [HI, which permits to obtain an analytical solution for the pair correlation functions. The 
MSA correlation functions can be used to evaluate an equation of state [l, 0] • ln spite of quite 
reasonable results obtained for the structure and thermodynamics in a general sense, neither virial 
nor compressibility equations of state permit to obtain gas-liquid transition for the RPM, which 



contradicts the computer simulations [8Hl5j. By contrast, the equation of state resulting from the 
energy route predicts the first-order phase transition between an ionic vapor and a dense ionic 
liquid @, 0|- However, there is a big discrepancy between the values of critical temperature and 
critical density of MSA theory and the critical parameters resulting from computer simulations. 
Moreover, a more sophisticated generalized mean-spherical approximation [l6j], which was designed 
to reconcile thermodynamic inconsistency of the MSA, does not lead to a better agreement of the 
critical parameters with simulations as well. 

Since both MSA and generalized MSA theories underestimate the critical density and overes- 
timate the critical temperature, attempts were made to include ad hoc the effect of ion pairing 
(association) into the RPM, which is most pronounced along the vapor branch of the coexistence 
envelope and near the critical point [l^, [13] ■ The inclusion of the association into the theory 
improves the agreement of theoretical predictions for the critical constants with simulations (cf. 
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table I of [13]). However, the overall shape of the coexistence envelope coming from several theoreti- 
cal approaches is not well described compared to simulation data. In particular, the liquid densities 
along coexistence are poorly described for different temperatures. We should bear in mind that 
the Hamiltonian of the model with association used in the theoretical developments differs from 
the genuine RPM used in simulations. Recent developments in the modelling of different types 
of associations and clustering in liquids and solutions have been comprehensively and critically 
reviewed by Holovko in |18l ] . 

The liquid-vapor interface is one of the simplest examples of nonuniform systems. However, 
compared to the case of Lennard- Jones fluids [19], much less is known about the liquid- vapor 
interface in the RPM. Theoretical work on the gas-liquid interface for RPM was pioneered by Telo 
da Gama et al. [2(|, that used the gradient expansion method [2lj]. The problem of describing 
an ionic liquid- ionic vapor interface has been also studied by Groh et al. [22], who proposed a 
density functional approach to evaluate the surface tension and the density profiles through the 
interface. Their approach involved a local density approximation for the hard-sphere part of the 
free-energy functional and a nonlocal treatment of Coulombic contributions. The latter term was 
evaluated approximating the inhomogeneous pair correlation functions by their bulk counterparts. 
There are significant differences between the theory used by Groh et al. and the density 
functional approaches that were proved successful in the description of electrical double layers at 
charged walls [23M3lj | . Namely, the theories of (UH^I treat the electrostatic part of the free energy 
by means of a second-order density expansion about the density of a reference fluid, which was 
taken to be the homogeneous bulk fluid far from the surface. Although such an approach can be 
appropriate for many purposes, but it is problematical when it comes to a liquid-vapor interface, 
where two bulk fluids are involved. Moreover, it is known that in the case of Lennard- Jones fluids, 
the corresponding second-order expansion with respect to a homogeneous fluid fails to account 
for the liquid-vapor coexistence. More recently, the description of a liquid-vapor interface was 
the subject of investigations of Weiss and Schroer [32| . Using a square-gradient type theory they 
computed the density profiles and the interfacial tension at different temperatures using Debye- 
Huckel theory and its extension to ion-pair formation, as well as adding the free energy term 
descri bing correlations between ion pairs as entities and free ions, as developed by Fisher and 
Levin [17j . 

Any density functional theory determines thermodynamic properties of an inhomogeneous fluid 
from the Helmholtz free energy F and its functional dependence on the local densities {pi(r)}. The 
free energy functional is commonly decomposed into the sum of three contributions, namely into 
the ideal, the hard-sphere, and the electrostatic terms. Of course, one can also include here an 
additional term, due to possible association of ions. Various formulations of the DFT have been 
discussed in the literature. As we have already stressed, a majority of the DFTs have followed 
perturbational second-order exp ansion of the electrostatic free energy functional with respect to a 
bulk homogeneous fluid [23" 



Recently, Gillespie et al. 33|, 34| proposed a version of the electrostatic free energy functional 
that replaces a uniform reference system with a suitably chosen position-dependent reference fluid. 
The inhomogeneous reference fluid densities are then computed from the local densities by a self- 
consistent iteration procedure. Actually, Gillespie et al. [H, [34J proposed a reference fluid density 
functional which permits to construct a reference model that locally satisfies electroneutrality con- 
dition and has the same ionic strength at every point as the inhomogeneous fluid in question. Such 
a construction of the reference fluid permits to apply the expression for the electrostatic contribu- 
tion to the free energy, which results from the equation of state for a bulk ionic system and makes 
it unnecessary to employ the direct correlation functions as input quantities. This kind of approach 
was proposed by us to study liquid- vapor transitions in RPMs confined to slit-like pores [3ll [35j . 
The results of our approach reasonably well reproduce the structure of ionic fluids at a charged 
wall at sufficiently high temperatures and correctly predict (3(| the temperature dependence of 
the double layer capacitance, in agreement with experiments and simulation data [37| . Quite re- 
cently, the theory was also extended to include the effects of association between unlike ions [Hj]. 
However, the applications of the approach described above would not be complete without explor- 
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ing liquid-vapor interface of ionic fluid. Therefore, in this communication we intend to apply the 
theory of |35l . l36l [38j to the study of liquid-vapor interface of the RPM. We calculate the density 
profiles and the interfacial tension. Moreover, we would like to investigate the effect of association 
leading to the formation of ion pairs and whether this type of effects improves the description of 
the interfacial properties similarly to the bulk structure and thermodynamics. 



2. Theory 

We consider a binary mixture of ionic species. The symbols d, Zi = qi/e and pi denote, respec- 
tively, the hard-sphere diameter, valence of ions and the chemical potential of species i = 1,2, and 
e is the electron charge. The interaction between the ions is 

n _ / oo, r < 1, , , 

Uij[r)-<y JZiZj/er, r > 1, (ij 

where r is the distance between a pair of ions. We also assume that Z\ = — Z 2 = 1 and that the 
dielectric constant e is uniform throughout the entire system. 

The DFT we use in this work is identical to that described in [35 , 38 1 and therefore we present 



only its brief description. If there is no external potential field the grand potential of the system is 
written in the form, 

n= J drf[{pi}]- J2 I Pi(z)»i&- (2) 

J 1=1,2'' 

According to the usual density functional treatment, the free energy density functional, /[{pi}], 
is decomposed into ideal, hard sphere, electrostatic and associative terms /[{/Ci}] = /id[{pi}] + 
f hs[{ Pi}] + fd[{pi}] + faa[{Pi}]- The ideal term is fid[{pi}] = X) j= i )2 [PiO*) ln Pi( z ) ~ M Z ) L w hereas 
for the hard-sphere contribution w e ap ply an expression resulting from a recent version [39( of the 
Fundamental Measure Theory |4oL Ulj . with the free energy density consisting of terms dependent 
on scalar and on vector weighted densities, n Q (r) (a = 0, 1, 2, 3) and n Q (r) (a — VI, V2) [39j], 

, W1 s , nm 2 -n vl -n V2 3 , 2 ,n 3 + (1 - n 3 ) 2 ln(l - n 3 ) 

/(hs) = -no In 1 - n 3 ) + + n$(l - 3^ — — w ^ , 3 

v ' 1 — n 3 dbnn^l — n 3 ) z 

where ^(r) = |ny2( r )|/"2(r). The definitions of weighted densities, n Q (r), a = 0, 1, 2, 3, VI, V2, 



are given in [40|, |4l| . 



The electrostatic contribution is evaluated from the approach described in [35, 38 1, according 
to which 



e 



2 



nipt}) , r 3 ({p 4 }) 



umvxr = + ^ TT f(md + ^ ' (4) 

where pi(z) denote suitably defined inhomogeneous average densities of a reference fluid. The form 
for fe\\\Pi(z)}] results from the MSA equation of state evaluated via the energy route 0, 0, [B~ 
I n the above ^({pi}) — + 2K({pi})d — l)/2d, where n({pi}) denotes the inverse Debye 
screening length, K 2 ({p,}) = {Aire 2 jT*)\Z\p x + Z 2 p 2 ]. 

The reference fluid densities {pi} are evaluated by employing the approach of Gillespie et 
al. [33L [33|. In the case of gas- liquid interface, the spatial symmetry of the RPM implies that the 
density profiles of the two ionic species should be the same. Therefore, in the presently considered 
situation, the approach of Gillespie et al. [33, 34 1 is simplified and the reference fluid densities are 



given by 

Pi(z)= I Pi (z)Wi(\r-r'\)dr', (5) 



where the weight function Wi(r) is just a normalized step function 



fi (| r ^|) = f^W) ( 6 ) 
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The radius of the sphere over which the averaging is performed, i?/, is approximated by the 
"capacitance" radius, i.e., by the ion radius plus the screening length 



Rf(r) = - 



1 



2r(te(r)}) 



(7) 



The evaluation of Rf requires an iteration procedure. This iteration loop should be carried out in 
addition to the main iteration procedure used to evaluate the density profiles, see for details [33l l3~il | . 

Finally, the associative contribution, / as , is formulated at the level of the first-order thermody- 
namic perturbation theory |38l . l42j 



fas/kT = ^2 Pi 
i=1.2 



,r 1 Ot({pi\) 

\na({p t }) + --^^-L 



(8) 



where a is the degree of dissociation according to the mass action law, 2a = 1 — K(pi + p-i)® 2 ■ 
The association constant K = K({pi},T) is a product of solely temperature dependent term, K° 
and iT?, K = K°K^. The constant K°, is 



8nd 3 Y, 



(T* 



-2m 



^ (2m)!(2m - 3) ' 



m=2 



whereas K 1 follows from the so-called simple interpolation scheme [l3, H, 44 1 



K , = a - m 



(1 - fj)3 



exp 



Td {2 + Td) 
T* (1 + Td) 2 



(9) 



(10) 



The theory presented above uses different weighted densities to evaluate the hard sphere and 
electrostatic free energy contributions. Note that the idea of using different weighted densities to 
evaluate these contributions has also been employed by Patra et al. [2R [26| . 

The density profile is obtained by minimizing the excess grand potential functional AQ — fl—flb, 



SAfl 

Spi(z) 



0, for = 1,2, 



(11) 



where fij, is the grand potential of a bulk uniform system of density {pb,i}- 

The evaluation of the density profiles requires the knowledge of the densities of both coexisting 
phases. Therefore, we have precisely evaluated the bulk phase diagram prior to the density profile 
calculations. Next, the local density calculations were carried out assuming that for z ^ — L z 
and for ^ L z the density of the fluid is equal to the bulk densities of gaseous and liquid phases, 
respectively. The value of L z was evaluated at each temperature by testing the convergence of the 
density profiles towards the final solutions. We have started with L z = 30d and the consecutive runs 
were carried out doubling the value of L z for each new run. All the integrations were performed 
using Simpson method with the grid size of O.Old. The convergence criterion definitly states that 
the maximum difference between two consecutive iterations must be smaller than 1E-8 percent. 
The knowledge of the density profiles of ions permits to calculate the excess grand thermodynamics 
potential per unit surface area, A, i.e., the surface tension 7 = Ail/A. 



3. Results and discussion 

Figure 1 shows the phase diagram in the density-temperature plane. The phase diagram result- 
ing from the MSA energy equation of state has been presented in several previous works, cf. fl3ll2^. 
Nevertheless, we have decided to include it here for completeness of the study. The solid line corre- 
sponds to the system without association (i.e. the associative free energy term is neglected in the 
grand thermodynamics potential functional), whereas the dashed line describes phase behavior of 
the model with the association effects included. Inclusion of chemical association into the theory 
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provides a mechanism for the formation of ionic pairs. The fraction of pairs depends on density, on 
temperature and on the association constant. From physical point of view, the formation of pairs 
of oppositely charged ions, alters the effective interactions between all the particles. In effect, the 
critical temperature of liquid-vapor transition, T* « 0.0745, becomes lower compared to the criti- 
cal temperature of the model without association, T* s» 0.0786. The liquid-vapor phase diagrams 
are strongly asymmetric. At low temperatures the vapor phase is very dilute. For example, for the 
system without association at T* = 0.05 the vapor density is p'l w 1.094E — 5, whereas the density 
of the coexisting liquid phase is pi « 0.2135. 




Figure 1. Phase diagram resulting from MSA approach and the energy equation of state. Solid 
and dashed lines denote the results for the model without and with association. 




Figure 2. A comparison of the total density profiles, p* (z) — (pi (z) + pi_(z))d 3 resulting from the 
DFT without association theory (solid lines) and from MSA1 theory of Groh et al. [22] (dashed 
lines). The temperatures T/T c are given in the figure. 



First, we consider the system without association. The evaluation of the liquid- vapor density 
profiles was carried out assuming that the limiting values pi(—L z ) and Pi(L z ) are equal to the MSA 
densities of coexisting phases. Of course, the agreement of our results with the results published 
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previously [22j, |32| essentially depends on the agreement between the bulk phase diagrams. One 
should keep in mind the above remark while analyzing figure 2, where we show a comparison of our 
results with those of Groh et al.[22||. The latter profiles were obtained from a different bulk theory 
(called by them "the MSAl approach" [22]), which yields the value of the critical temperature, 
T* « 0.0846, different from that obtained by MSA theory. Note that a comparison of the entire 
phase diagrams resulting from the MSA and MSAl approaches in given in figure 1 of [22j. 

; — 1 — i — 1 — i — 1 — i — : 

• •*:>-,, T*= 0.074 




2 4 

z/d 



Figure 3. The convergence of the liquid parts of the single species density profiles towards 
the bulk liquid equilibrium density at three different temperatures given in the figure. Figure 
also contains the equations for straight lines (marked by solid lines) approximating the data, 
computed from density functional theory. Here y = \n\p)b, i 1 — pi(z)] and p\ A is the single species 
density of an ionic liquid at coexistence. The calculations are for the model without association. 



The way in which the density profile approaches the asymptotic bulk values has been discussed 
by Groh et al. [22j. We have checked that the decay of the function p\ i — p%{£) (where p\ i is the 
liquid-phase bulk density) with z — > L z is exponential, i.e. ln[pi(z) — p b J oc z/£. This point is 
illustrated in figure 3, where we have plotted ln[p b i — Pi(z)] versus z/d. At larger distances this 
dependence is fairly well approximated by a straight line with the correlation coefficient, R 2 , very 
close to unity. 

We now consider the effect of the association, cf. figure 4. Part a shows a comparison of the 
total density profiles for the systems with (symbols) and without association (lines) at different 
temperatures. At low temperatures the liquid density for the system with the association included 
is higher than for the system with the association effects switched off, cf. figure 1. The profiles for 
the model with association are more "smeared out", but this is the effect of lowering the critical 
temperature due to association effects. Similarly to the case of a system without association we 
have found exponential decay of the function p\ i — pi{z) with z — > L z (the corresponding plot 
was omitted for the sake of brevity). Figure 4b shows the profiles of a(z) (the phase diagram in 
the a-T plane has been already presented in [38j]) We stress that the interfacial region of a(z) is 
shifted towards rarefied phase compared to the density profiles. It means that the position of the 
"pseudo-Gibbs" dividing surface, which can be introduced for the dissociation degree in analogy 
to the usual liquid-vapor dividing surface is also shifted compared to the position of the usual 
dividing surface (which is located at z = 0) . Note that similar shift of the dissociation degree was 
also observed in the case of liquid- vapor interfaces of non- ionic fluids [45| . 

The surface tension obtained from the density functional theory and from the MSAl and MSA 
approaches, developed by Groh et al. [22j], as well as from the square gradient theory [2(J are shown 
in figure 5. We have plotted here the ratio 7c? 2 /kT c versus the temperature reduced by the critical 
temperature, T/T c . Except for the temperatures very close to the critical temperature, the present 
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0.4 0.6 0.8 1 

T/T 

c 



Figure 5. The dependence of the surface tension on the temperature. Curves labelled MSA, 
MSA1 were evaluated by Groh et al. the curve labelled "square gradient theory" was eval- 
uated by Telo da Gama et al. fiol ]. the curve abbreviated as DFT was obtained for the model 
without association, whereas black circles - for the model with the association effects included. 



theory predicts lower values of the surface tension than all the rest approaches. The associative free 
energy terms still lower the value of 'yd 2 /kT c . Note that the MSA theory grossly underestimates 
the densities of the coexisting liquid phase. This can suggest that the values of the surface tension 
are also underestimated. On the other hand, the value of the critical temperature is overestimated 
by the MSA approach, and hence can lead to too high values of the interfacial tension. 

Near critical temperature, both models (i.e., the model with and without association) yield 
the standard mean-field critical behavior, i.e., ^yd 2 /kT c oc (1 — T/T c )( 3 / 2 ' (see figure 6). As the 
temperature approaches the critical point, the theory yields a nearly linear variation of In [7 d 2 / kT c ] 
with ln[(l — T/T c )\ and the slope of the straight line approximating the obtained data is almost 
perfectly equal to (3/2). 
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f Yd 2 /kT c oc(l-T/T c ) (3/2) 



in t I I I I I I I I I I L. 

-6 -5 -4 -3 -2 -1 
ln(l-TYT) 



Figure 6. The decay of the surface tension on approaching the critical temperature. Solid line 
and empty circles are for the model without association, whereas dashed line and filled circles - 
for the model with association. The straight lines slope is 3/2. 



3 - 



T3 

l 2 

a 



n 1 1 1 1 1 1 1 1 1 r 

W/doc (l-T/T c ) ( " 5) 




-6 



° • 

o « 



_L 



_L 



-4 -3 -2 
In(l-TTT) 



Figure 7. The divergence of the interfacial width W , on approaching the critical temperature. 
The straight lines slope is —1/2. Abbreviations are the same as in figure 6. 



The width of the interfacial zone can be characterized by the parameter W, denned as [46 



W 



-\p(z = L z ) - p(z = -L z )\ 



dp(z) 



dz 



(12) 



where zq is given by p(zq) = (l/2)[p s (z — L z ) + p(z = — L z )] and p(z — L z ) and p[z — —L z ) are the 
total densities of the coexisting liquid and vapor phases. Obviously, the interface becomes wider 
as the temperature increases. The effects of association slightly increases the interface width. As 
the temperature approaches the critical temperature, the interfacial width diverges. To investigate 
the character of this divergence, we have plotted the values of \nW/d versus ln(l — T/T c ), see 
figure 7. For T approaching the critical temperature, the dependence of In W/d versus ln(l — T/T c ) 
is linear. We have found the following scaling: W* oc (1 — T/T C (M))~ 5 . The obtained value of 
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the exponent (0.5) is characteristic of mean-field type theories [22j. The scaling of W is consistent 
with the scaling of the surface tension. 




Figure 8. A comparison of the total density profiles for the models without association (solid 
and dotted lines) and with association (dashed and dash-dotted lines) with computer simulation 
data. The simulation results for the profiles (see details in the text) are shown by symbols [47ll4sl |. 
The temperatures are given in the figure. 



Our discussion involved the results of different, but solely theoretical approaches. However, it 
must be clarified how accurate our approaches are (with and without association) compared to the 
available computer simulation data. In particular, Alejandre and co-workers investigated the liquid- 
vapor interface for the RPM using molecular dynamics simulation method supplemented by certain 
specific technical peculiarities dealing specifically with the RPM (47, 4S|. These authors presented 
the density profiles, p*(z) at four reduced temperatures, T*=0.035, 0.038, 0.040 and 0.043 lijl . On 
the other hand, a single density profile at T* — 0.391 relevant to our study was given in [48] (this 
profile is for the soft primitive model but with parameters permitting comparison with the RPM). 
The authors have not performed estimates of the critical parameters coming from their method 
of study of the RPM. Thus, they somewhat arbitrarily used the critical temperature obtained 
in T* = 0.0492, for their purposes. We have rescaled the temperature and density with 
respect to the critical parameters of both theoretical approaches (with and without association) , 
respectively, and present here a comparison of the density profiles with computer simulation data 
in figure 8. Note that the critical density given in was used to rescale the simulation data. 
We observe that the rescaled liquid phase density at coexistence coming from the DFT for the 
model with ion pairing is in reasonable agreement with the simulation data. The DFT for the 
model without association greatly overestimates these liquid densities. However, the width of the 
intcrfacial region is not described well. Theoretical approaches yield narrower interface width for 
two reduced temperatures in question. It is difficult to attribute this specific inaccuracy to the 
particular term of the free energy functional. However, it seems that the mean field consideration 
requires improvement in order to better describe interparticle correlations and hopefully to obtain 
a more adequate interface width. 

As concerns the relationship between the simulation data for the surface tension and the results 
of theoretical approaches we would like to comment on the following points. A set of computer 
simulation results for the surface tension were put together with the predictions of some theoretical 
approaches [2(| 22, 3^1 (that were discussed in the introductory part of our communication) on 



arbitrary temperature scale, see figure 2 of [47|. Apparently, the theories yield the surface tension 
of the same magnitude as computer simulation data. This permitted to conclude that the surface 
tension from this set of theories 2^, 22 1 is overestimated because the effect of ion association is not 



13603-9 



A. Patrykiejew, S. Sokotowski, O. Pizio 



taken into account and solely the theory that accounts for association [32] is in good agreement 
with simulations. Then, it would seem that our approach taking ion pairing into account is the 
best, see figure 5. Moreover, a comparison of the theories, simulations and experimental data was 
performed [48] using an arbitrary scale. Such an idyllic picture is, however, entirely destroyed, if one 
rescales the reduced temperature scale by the critical temperature of each theory and of simulations. 
This is actually an adequate comparison. According to [4j| the simulations yield, e.g. 'yd 2 /kT c =0.3 
at T/T c =0.6. Such a value is very much higher than the results of any theory, c.f. figure 5, at 
this particular temperature. Moreover, this means that the inclusion of association effects makes 
things worse than from any other theory. This is difficult to accept, however. According to the 
more recent simulation work of the same authors (48| it is difficult to draw a definite conclusion 
about the dependence of the surface tension on temperature for the model in question due to the 
statistical error associated with the value for surface tension. Moreover, it has been discussed [48| 
that the inflection point at T* — 0.04 (again it is a reduced temperature of simulation) can be 
present on the dependence of surface tension on temperature as a result of chemical association 
of ions in the RPM in close similarity to hydrogen bonded fluids. In our opinion, additional work, 
both theoretical and simulational, is necessary to obtain definite answers about the behavior of 
surface tension on temperature for the model in question and the related models. From theoretical 
viewpoint one needs a theory that yields a more adequate shape of the bulk coexistence envelope. 
Then, the free energy expression coming from this approach, if available, would contribute to the 
development of better density functional approaches for inhomogeneous ionic fluids. 

To summarize briefly, in this work we have applied density functional theory to the study of 
the liquid-vapor interface of a RPM fluid. Of course, the theory we use here is ad hoc. However, 
our previous calculations (35l . l36| have shown that it reproduces reasonably well the structure 
of a fluid at a charged and uncharged wall, predicts the dependence of the capacitance of the 
double layer on the temperature and yields "capillary evaporation" phase diagrams for confined 
ionic systems [35, 36|, 38|. We have shown that the inclusion of ion pairing effect leads to a better 
agreement of the density profiles of ions across the liquid-vapor interface. A more sophisticated 
approaches (l3| . compared to the present one, including the effects of ionic association can lead 
to even better agreement of the coexistence envelope with simulations and of the density profiles 
across the liquid-vapor interface. We plan to extend our study in this respect in the nearest future. 
However, the inaccuracies of the values for the liquid phase density at coexistence at different 
temeperatures prevent us from obtaining a reasonable agreement with the available simulation data 
for the surface tension. It seems, however, that a more ample set of simulation data using different 
techniques is necessary. Nevertheless, the theories of this study yield a decreasing surface tension 
with increasing temperature as one intuitively would expect. Still we are not able to conclude that 
the theoretical approach is entirely successful both for the fluid-solid and fluid-fluid interfaces and 
look forward to improving it in future work. 
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Mi>Kcpa3Ha rpaHMi45i pi^i/ma-napa o6ivie>KeHoV npi/iMiTMBHoY MOfle/ii 
ioHHMX n/iMHiB 3 TeopiT cpyHKi4ioHa/iy rycTi/mn 

A. riaTpnKee^, C. Coko/iobckP, 0. l1i3icP 

yHiBepcuTeT iM. MapiT CmiaflOBCbKoT-Kiopi, JlKi6jiiH, Pecny6/iiKa l~lo.nbLU,a, 
iHCTMTyT xiiviiT yTAM, KowoKaH, MeKCWKa 

M 1/1 flocnifl>KyeMO Mix<cpa3Hy rpa hm l^kd piflMHa-napa o6Me>KeHoT npi/iMiTMBHoT MOfle/ii ioHHoro n/iwHy, bhko- 
pucTOByKDHH Teopiio cpyHKLiioHajiy rycTi/iHi/i. 3acTocoBaHa Teopis BK/nonae e/ieKTpocTaTM h hm m BK/iafl y cpyH- 
Ki4i0Ha.fi Bi/ibHoT eHepriT, jrkmW Bi/iHi/iKae 3 eHepreTi/iHHoro piBHflHHfl CTaHy fl/ia o6'eMy, BKJiafl Bifl cepeflHbo- 
ccpepuHHoro Ha6jiw>KeHH5i fl/isi o6Me>KeHoT npi/iMiTHBHoT MOfle/ii, a Ta ko>k accmiaTWBHwti BK/iafl, skhB BMHWKae 
b pe3y/ibTaTi BpaxyBaHHa yTBopeHH^ ioHHi/ix nap. Mm nopiBHioeMO npocpi/ii rycTHHH i 3HaneHHa noBepxHeBoro 
Ha~rary 3 pe3y/ibTaTaMH nonepeflHix TeopeTi/iHHi/ix niflxofliB. 

KyiKJHOBi enema: cpyHKU,ioHan rycTtAHt/i, ap,cop6u,\a, j}3hu,io>kkia, KpwcTanw 
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